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ABSTRACT 


Numerical  results  of  the  Ostwald  Ripening  problem  with  coalescence  are 
reported  for  aerosol  with  the  continuum  diffusive  growth  laws.  It  was  shown 
that  N (t)  -  t  \  where  N(t)  is  the  total  number  concentration.  These  re¬ 
sults,  though  consistent  with  the  asymptotic  theory  of  Ostwald  Ripening 
according  to  Li f shi tz-Slyozov  and  with  the  similarity  theory  for  the  continuum 
Brownian  coalescence  process,  show  that  the  scaling  arguments  are  insufficient 
to  establish  which  of  the  two  mechanisms  is  dominant  asymptotically.. 
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'Numerous  studies  have  been  concerned  with  the  asymptotic  stages  of 

growth  by  condensation/evaporation  of  droplets  or  particles  following  their 
appearance  by  a  homogeneous  nucleation  process.  Comparatively  less  atten¬ 
tion  has  been  given  the  associated  problem  of  growth  when  both  condensation/ 
evaporation  and  coalescence  occur.  For  aerosols  we  present  the  first  quan¬ 
titative  study  of  the  Ostwald  ripening  problem  with  a  realistic  coalescence 
process. 

Me  have  investigated  the  isothermal  spatially  homogeneous  growth  of 
particles  suspended  in  supersaturated  monomer  vapor  and  inert  host  gas.  The 
process  is  described  by  the  evolution  equation  for  the  singlet  density  func¬ 
tion  n(x,t): 

Vi  ZXo  3  x/2 

+  __  [v(x>s)n(x,t ) ]  *  /  *  b(x-x‘  ,x')n(x-x'  ,t)n(x'  ,t)dx' 

x 


-  n(x,t)  /  .  b(x' ,x)n(x‘ ,t)dx' 


where  n(x,t)dx  is  the  number  of  particles  having  masses  in  the  range  x,dx  at 
time,  t.  y(x,s)  is  the  growth  law  for  a  particle  by  the  condensation/evapora¬ 
tion  process  and  b(x' ,x)  is  the  coalescence  rate  coefficient  for  two  particles 

* 

of  masses  x1  and  x.  x  is  the  mass  of  the  smallest  particle  in  the  popula¬ 
tion  and  s  is  the  supersaturation  ratio  of  monomer  vapor.  We  do  not  consider 
here  competitive  growth  effects,  ^  as  the  theory  for  this  is  lacking  when 
coalescence  is  significant. 


Equation  (1)  is  coupled  to  the  mass  conservation  equation  for  monomer 


vapor : 


Cy  ||  -  -  <i»(x,s)n(x,t)dx  -  x%(x*,s)n(x*,t) 
x 


where  Cy  is  the  equilibrium  bulk  vapor  concentration. 

Equations  (1)  and  (2)  obey  the  mass  conservation  law: 


*  (cvs  ♦  V  -  0 

where  »  /"A  xn(x,t)dx  is  the  mass  concentration  of  particles.  An  impor- 
x 

tant  additional  moment  of  n(x,t)  is  the  total  .number  concentration,  N(t)  * 

/°°A  n(x,t)dx.  It  follows  from  Equation  (1)  that: 
x 

44  =  Y(x*,s)n(x*,t)  “/"*/"*  b(x‘  ,x)n(x*  ,t)n(x,t)dxdx' 

at  XX 

We  have  studied  the  evolution  of  n(x,t)  by  numerical  simulation  using 
techniques  described  in  detail  elsewhere.  These  simulations  use  the  con¬ 
tinuum  diffusive  growth  laws: 


¥>(x,s)  ■  *»7rRDCv(s  -  exp(K/R)) 


where  x  *  (4ir/3)R3d,  0  is  the  binary  diffusion  coefficient  of  monomer  in  host 
gas  and  K  is  the  Kelvin  coefficient  for  the  dependence  of  vapor  pressure  on 
droplet  curvature,  d  is  the  mass  density  of  a  particle  with  radius  R.  The 
rate  coefficient  for  coalescence  is  that  appropriate  to  the  continuum 
Brownian  process: 


b(x*,x)  -  (2kT/3u)(x1/3  +  x'1/3)(-j73  +  — {73) 


where  kT  is  the  thermal  energy  and  p  is  the  viscosity  coefficient  of  the 
host  gas. 

Equation  (1)  may  be  put  into  nondimensional  form  through  the  substitu¬ 


tions:  p  *  x/x  ,  t  -  4ttDC  (3/4*dx  2)'^3t,  $  »  x  n(x,t)/N  .  x  is  the  initial 
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mean  particle  mass  and  Nq  the  initial  number  concentration  of  particles.  In 


these  variables,  Equation  (1)  becomes: 


^  [p1/3(s  -  exp(K'/p1/3))$(p,T)] 


a  /P/2  bf  (p-p1  ,p'  )$(p-p*  .xJ^Cp'.Tjdp*  -  #(p,t)  I  bj  (p' ,p)^(p' ,t)dp'  (3) 
P  P 


S~ 


whe  re 


kTNoxo2/3/3nyOCv  (3/***d ) 1  /3 


The  coefficient  a  is  a  measure  of  the  relative  rates  of  coalescence  to  conden* 
dat ion/evaporat ion.  In  Equation  (3),  for  o  -*■  »,  the  coalescence  process  be¬ 
comes  dominant  and  for  a  ■*  0,  only  the  condensation/evaporation  process 

3  12 

remains  with  results  previously  described.  ’ 

8  1 6 

A  question  addressed  by  others  ’  using  scaling  arguments  is  which 
mechanism  -  condensation/evaporation  or  coalescence  -will  be  dominant  for 
intermediate  values  of  a.  The  conclusion  from  these  analyses  is  that 
asymptotically  coalescence  will  dominate  at  higher  particle  mass  concentra¬ 
tions  -  that  is  for  larger  values  of  a.  We  show  here  that  this  is  not  the 
case  and  that  apparently  coalescence  can  only  be  dominant  asymptotically 
for  the  growth  laws  studied  for  vanishingly  small  rates  of  the  condensation/ 
evaporation  process.  This  is  illustrated  by  presenting  a  result  from 
numerical  solution  of  Equations  (1)  and  (2)  for  the  stated  growth  laws. 
Parameters  used  yield:  a  ■  0.31,  t  *  5-3  x  103t,  s(t“0)  *  10,  K/Rq  m  1.0. 

Figure  (1)  shows  the  evolution  in  time  of  the  moments  N,  of  the 
distribution  together  with  variation  with  time  of  the  interfacial  diameter. 

Op  ;  Dp  ■  2K/lnS.  After  an  initial  induction  period  ’  i  t  can  be  noted  that 
N (t)  -  t  This  dependence  is  the  same  predicted  for  asymptotic  continuum 
diffusive  growth  by  the  condensation/evaporation  process  according  to  the 
Lifshitz-Slyozov  (LS)  theory3  It  is  also  the  same  dependence  found  both 
numerically  and  by  similarity  theory^ for  the  continuum  Brownian  coalescence 
process.  The  increase  of  Op  with  time  follows  the  law:  Op  -  t  ,  also 
the  same  predicted  by  the  LS  theory  for  only  the  condensation/evaporation 
process.  Interestingly,  even  substantial  coalescence  rates  do  not  alter 
significantly  this  result  from  condensation/evaporation  similarity  theory. 


"'/'/V.'/V.'  .•  v  v  '  .• 


Since,  separately,  condensation/evaporation  and  coalescence  both  give  the 
result,  N(t)  -  t  for  the  continuum  di ffusive  cases,  scaling  arguments  are 
insufficient  to  establish  which  of  the  two  mechanisms  is  dominant  asymptoti¬ 
cal  ly. 

Figure  (2)  shows  the  contribution  to  the  total  rate  of  change  of  N, 
dN/dt,  by  each  of  the  two  processes  -  condensation/evaporation  and  coales¬ 
cence.  The  continuous  curve  is  that  for  the  ratio  of  coalescence  rate  of 
change  of  N  to  the  total  rate  and  the  dashed  curve  that  for  the  condensa¬ 
tion/evaporation  process  according  to  Equation  (1).  Clearly,  the  contri¬ 
bution  of  coalescence  to  the  change  of  N  decreases  with  time  while  that  for 
condensation/evaporation  increases.  Therefore  previous  conjectures  on  the 
dominance  of  either  of  the  two  processes  appear  to  be  incorrect  for  these 
continuum  diffusive  growth  processes. 
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FIGURE  CAPTIONS 


FIGURE  1.  Evolution  of  total  particle  number  concentration,  N,  and 
total  particle  mass  concentration  for  continuum  diffusive 
growth  with  simultaneous  processes  of  condensation/ 
evaporation  and  coalescence.  0^  is  particle  diameter  at 
condensation/evaporation  interface  of  particle  size  distri¬ 
bution. 

FIGURE  2.  Relative  contributions  of  coalescence  (solid  curve)  and 
condensat ion/evaporation  (dashed  curve)  to  rate  of  total 
time  change  of  particle  number  concentration,  N. 


